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ABSTRACT 


Nuclear DNA ITS and ETS sequences of 71 representatives from nine genera and 11 sections of the core Laureae were 
combined with a matrix of morphological characters, analyzed using maximum parsimony with both equally and successively 
weighted characters, and analyzed for Bayesian inference, minimum evolution by neighbor joining. and maximum likelihood 
inference for molecular data alone. The large genera Actinodaphne Nees, Lindera Thunb., and Litsea Lam. were polyphyletic, 
as were Lindera sect. Aperula (Blume) Benth. and Litsea sections Conodaphne (Blume) Benth. & Hook. f.. Cylicodaphne (Nees) 
Hook. f, and Tomingodaphne (Blume) Hook. f. In contrast, Veofitsea (Benth.) Merr. was monophyletic and terminal in a larger 
monophyletic lineage above an Actinodaphne grade. A major disparity exists between these molecular results and traditional 
morphology-based classifications within Lauraceae. These results suggest that the use of two- versus four-celled anthers for 
Laureae generic delimitation has resulted in polyphyletic or paraphyletic genera, and the character of dimerous versus 
trimerous flowers is of only limited phylogenetic value. Several of the major lineages in Laureae are supported by inflorescence 
morphology and ontogeny, with Laureae defined by short shoots with a vegetative terminal bud, splitting. into. thyrsoid 
(Actinodaphne and Neolitsea) versus racemose (Laurus, Litsea s. str., and Lindera s. str. and Lindera sect. Aperula), although 
there appear to be at least two different pathways to form the Laureae pseudo-umbel. Similarly, imbricate, early deciduous 
inflorescence basal involucral bracts defined an Actinodaphne—Neolitsea Parasassafras -Sinosassafrus clade. although within 
il, Neolitsea was defined in part by decussate, persistent bracts. Accordingly, our study indicates the need for caution in the 
use of morphology for assessing affinity in Laureae, as virtually all traditional morphological characteristics show high levels of 
homoplasy and/or reversal, but future research may help to resolve whether this indicates problems of homology or ontogenetic 
convergence, 
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The core Laureae comprise eight genera with 
approximately 500 species, mainly from tropical and 
subtropical Asia, making this region a center of 
generic and infrageneric diversification and distribu- 
tion for the Laureae (Li, 1995). The group is defined 
morphologically by its dioecious breeding system, 
mostly pseudo-umbellate inflorescences subtended by 
involucral bracts, and introrse anther cells in the third 
whorl (Li & Christophel, 2000; Li et al, 2004). 
However, a large portion of what we now regard as the 
Laureae, including Actinodaphne Nees, Dodecadenia 
Nees, Headaphne Blume, Laurus L., Lindera Thunb., 


Litsea Lam., Neolitsea (Benth.) Merr., Parasassafras D. 
G. Long, and Sinosassafras H. W. Li, were only placed 
together relatively recently (Rohwer, 1993; van der 
Werff & Richter, 1996). Although half of the genera in 
the Laureae are monotypic or oligotypic, Actino- 
daphne, Lindera, Litsea, and Neolitsea each possess 
over LOO species. 

Hooker (1890) and Li et al. (1984) divided Litsea 
into four sections, but these divisions were not 
explicitly phylogenetic. Section Litsea is evergreen 
with alternate, penninerved leaves, a racemiform 


inflorescence, and non-enlarged perianth tubes with 
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reduced or absent perianth lobes. Section Conodaphne 
(Blume) Benth. & Hook. f. is evergreen with alternate 
or opposite penninerved leaves and non-enlarged or 
slightly enlarged perianth tubes, whereas section 
Cylicodaphne (Nees) Hook. f. represents evergreen 
taxa with alternate, penninerved leaves, an enlarged 
perianth tube, and cup-shaped fruiting cupule. In 
contrast, section Tomingodaphne (Blume) Hook. f. has 
deciduous, alternate, penninerved leaves and non- 
enlarged, six-lobed perianth tubes. 

Tsui (1987) divided Lindera into eight sections, 
mainly following Hooker's (1890) criteria for Litsea, but 
also incorporating Li's (1985) concept of the "shortened 
brachyblast,” which had been proposed to explain 
inflorescence evolutionary trends in Laureae. Section 
Lindera has taxa with deciduous, penninerved leaves 
and well-developed terminal buds on shortened 
brachyblasts; section Sphaerocarpae H. P. Tsui pos- 
sesses deciduous, triplinerved or trinerved leaves and 
well-developed terminal buds on shortened brachy- 
blasts; and section Palminerviae Meisn. was created for 
deciduous taxa with lobed, trinerved leaves and well- 
developed terminal buds on shortened brachyblasts. 
Section Aperula (Blume) Benth. has evergreen, penni- 
nerved leaves, well-developed terminal buds on 
shortened brachyblasts, and long-pedunculate, racemi- 
form-arranged inflorescences; section Cupuliformes H. 
P. Tsui possesses evergreen, penninerved leaves, 
funnel-shaped glands on the third anther whorl, and 
enlarged perianth tubes forming cup-shaped fruiting 
cupules. Section Daphnidium (Nees) Benth. has 
evergreen, trinerved or triplinerved leaves, and non- 
developing terminal buds on shortened. brachyblasts; 
section. Polyadenia (Nees) Benth. is evergreen with 
penninerved leaves and ill-developed terminal buds on 
the axillary short shoots, whereas section Uniumbella- 
tae H. P. Tsui has evergreen, trinerved leaves, and the 
long-pedunculate pseudo-umbel is solitary and borne 
in the axil of a normally developed leaf. 

Generic delimitation in Laureae is problematic 
1957; Hutchinson, 1964; Richter, 


1981), the major point being the significance of two- 


(Kostermans, 


versus four-celled anthers as generic descriptors, 
although recent Lauraceae classifications tend to 
downplay this feature (e.g, Rohwer et al., 1991; 
Rohwer, 1993; van der Werff & Richter, 1996). 

The most recent hypothesis about evolution within 
the Laureae was constructed under the assumption of 
maximum parsimony (MP) using nucleotide data from 
chloroplast DNA (cpDNA) matK and nuclear DNA 
(nrDNA) ITS (Li et al., 2004). This study included 22 
species and one outgroup taxon, and shed light on 
major lineages, as well as on the relationships of some 
previously problematic genera. That study suggested 
that Litsea and Lindera are polyphyletic, with no 
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support for characters such as two- versus four-celled 
anthers, and identified novel clades that appeared to 
reflect inflorescence structure and ontogeny. Never- 
theless, although the molecular analyses provided 
some phylogenetic insights of the core Laureae, 
generic boundaries and relationships were still 
uncertain due to limited sampling and the fact that 
morphology was not included explicitly in the 
analyses. In addition, the level of matK divergence 
was particularly low in Laureae, resulting in poorly 
supported groupings. 

To reconstruct a more reliable phylogeny of the 
core Laureae, it was necessary to sample additional 
taxa, as well as to utilize more informative DNA 
regions. The most widely used sources of molecular 
data for plant phylogenetic studies at the specific and 
generic levels are the [TS-1 and ITS-2 regions in 
nrDNA (Baldwin et al., 1995). However, for Lauraceae 
(Chanderbali et al., 2001) and Laureae (Li et al., 
2004), 


uninformative to provide robust node inferences by 


these two spacers were too invariant or 
themselves. In contrast, although generally less easy 
to sequence than ITS, the nrDNA ETS has been shown 
to be useful for supplementing ITS data in plant 
phylogenetic studies (e.g., Baldwin & Markos, 1998; 
Clevinger & Panero, 2000; Markos & Baldwin, 2002; 
Plovanich & Panero, 2004), giving it the potential to 
provide additional characters in the Laureae. 

In addition, because most recent molecular studies 
in Lauraceae have not explicitly included morphology 
as part of the analysis, preferring instead to discuss 
character patterns post hoc in relation to the 
molecular trees, there is a need for a combined 
analysis that includes both data sets in order to assess 
the impact, if any, of morphological data on the 
results. 

Accordingly, we have expanded the ITS study of Li 
et al. (2004) to include 71 taxa from the Laureae, 
three outgroups, and partial sequences from the ETS 
region, as well as information from morphological 
characters considered to be useful or definitive above 
the species level within the Laureae. The main goals 
were to assess the robustness of previously defined 
clades that had shown weak bootstrap support, as well 
as to see if support for features such as inflorescence 
structure was maintained with more extensive sam- 
pling and followed the explicit inclusion of morphol- 
ogy into the analyses, rather than just post hoc by 
character mapping. In particular, the hypotheses 
being tested were (1) Laureae genera based on two- 
versus four-celled anthers are not phylogenetically 
supported; (2) Litsea, Lindera, and Actinodaphne are 
not monophyletic; (3) Neolitsea is monophyletic; and 
(4) inflorescence ontogeny is phylogenetically impor- 
tant in the Laureae. 
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METHODS 
PANON SAMPLING AND DNA EXTRACTION 


In addition to 21 samples used previously by Li et 
al. (2004), except for Lindera tienchuanensis W. P. 
Fang & H. S. Kung (sect. Üniumbellatae) and L. 
thomsonii €. K. Allen (sect. Daphnidium) for which 
there was no more available material. an additional 50 
species were included to cover the sections within 
Lindera and Litsea variously recognized by Li et al. 
(1984), Bentham (1880), Hooker (1890), Kostermans 
(1957), and Tsut (1987), as well as to sample more 
(Table 1). 


sampling included one species each from Cinnamo- 


extensively within Neolitsea Outgroup 
mum Schaeff., Sassafras J. Presl, and Umbellularia 
(Nees) Nutl., as these genera were sister to the core 
Laureae in Chanderbali et al. (2001) and Li et al. 
(2004). 

Total genomic DNA was isolated from silica gel- 
dried material or herbarium specimens following 
Doyle and Doyle (1987), as modified by Li et al. 
(2004). Because of suspected fungal contamination in 
the previous ITS sequences of Sinosassafras flaviner- 
via (C. K. Allen) H. W. Li and Neolitsea confertifolia 
(Hemsl.) Merr. (Li et al. 2004). both 


extracted. 


were re- 


AMPLIFICATION, SEQUENCING, AND SEQUENCE, ALIGNMENT 


The [TS and 5.08 regions were amplified following 
the general methodology of White et al. (1990) and 
polymerase chain reaction (PCR) primers of Chan- 
derbali et al. (2001), using the minor modifications 
reported in Li et al. (2004). The complete ETS 
intergenic spacer (IGS) was amplified with the 18S- 
IGS and 265-F primers following the protocols of 
Baldwin and Markos (1998) and Starr et al. (2003) 
using Actinodaphne trichocarpa C. K. Allen, Lindera 
latifolia Hook. f., Lindera megaphylla Wemsl.. Litsea 
cubeba (Lour.) Pers., and Neolitsea chuii Merr. This 
procedure yielded a product between 2 and 3 kb. of 
which ca. 650 bases were sequenced using the 185- 
IGS primer. A conserved region located about 450 
bases from the ETS/188 boundary was selected to 
design a primer, ETS-1 (CCA GAA CTC GCA CTT 
GCT GAG CTT), and this region was amplified for all 
laxa using ETS-I and 188-08 as primers. PCR 
amplification followed the short-distance PCR proto- 
col of Baldwin and Markos (1998), with an annealing 
temperature of 55 C. and TaKaRa ExTaq (TaKaRa 
Biotechnology, Dalian, China). ALE PCR amplifications 
included negative controls to detect contamination. 

Although most PCR amplifications resulted in a 
single band, some individuals produced two size 


classes of PCR product. To isolate each of the two 


bands. PCR products were purified and blunt-end 
ligated into the EcoRV sites of the pMD18-T Simple 
Vector (TaKaRa) using the Original TA Cloning Kit 
(available from TaKaRa), and the fragments obtained 
following the manufacturer's protocol. PCR. products 
were purified using the QlAquick. PCR purification 
Kit (Qiagen, Tokvo. Japan) 
provided by the manufacturer, and sequenced in both 


following protocols 


directions. Sequencing reactions used the same 
primers as amplification and were conducted on an 
Applied Biosystems (Foster City, California, U.S.A.) 
3100 DNA automated sequencer. 

Sequences were checked against GenBank for the 
taxa of closest matching sequences by means of the 
BLAST search. Sequence chromatogram output files 
were aligned initially and edited base by base with the 
program Zei Man H (Lasergene software package; 
DNASTAR Ine., Madison, Wisconsin, U.S. A.). The 
edited sequences were then realigned using ClustalX 
version. 1.8 (Thompson et al.. 1997), and additional 
manual adjustments were made, if necessary, using 
BioEdit (Hall, 1999), AIl sequences analyzed in this 
study were deposited in GenBank (Table 1). 


VORPHOLOGY 


Morphological data were collected from personal 
observation of live plants, herbarium specimens, and 
literature. Nineteen morphological 
(Table 2) 


evolution and mapping because they have been 


the available 
characters were chosen for character 
regarded as important in delimiting genera within 
the Litsea complex, and/or between sections within 
Lindera or Litsea (Bentham, 1880: Hooker, 1890; 
Kostermans, 1957: Li et al., 1984; Li, 1985; Tsui, 
1987: Rohwer, 1993: van der Werff & Richter, 1996; 
Li & Christophel, 2000: van der Werff, 2001: Li et al., 
2004). A 


matrix of 19 discrete morphological 


characters was constructed for all of the taxa 


investigated (Table 3). 


PHYLOGENETIC ANALYSES 


We used several methods to reconstruct phyloge- 
nies and assess branch support for clades: MP with 
hoth equally and successively weighted characters. 
bootstrapping. Bayesian inference, minimum evolu- 
(ME) by 


likelihood (ML) inference. In this way, clades with a 


tion neighbor joining, and maximum 
strong. consistent phylogenetic signal should be 
recovered regardless of the methodology. 

MP analyses were performed using PAUP* version 
4.01b 10 (Swofford, 1999). Congruence of the ITS and 
ETS data sets was assessed using the partition 


homogeneity test (PHT; Farris et al. 1995) as 
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Table l. 


List of samples, voucher collection information, and GenBank accession numbers. 


Taxon 


Actinodaphne Nees 


A. cupularis (Hemsl.) Gamble 


A. forrestii (C. K. Allen) Kosterm. 

A. henryi Gamble 

A. kweichowensis Yen C. Yang & 
P. H. Huang 

A. lecomtei C. K. Allen 

A. obovata (Nees) Blume 

A. omeiensis (H. Liu) C. K. Allen 


A. paotingensis Yen C. Yang & P. H. 


Huang 
A. pilosa (Lour.) Merr. 
A. trichocarpa C. K. Allen 
A. tsaii Hu 


Dodecadenia Nees 


D. grandiflora Nees 


Iteadaphne Blume 
JI. caudata (Nees) H. W. Li 


Laurus L. 


L. nobilis L. 


Lindera Thunb. 


Sect. Aperula (Blume) Benth. 
L. latifolia Hook. f. 
L. longipedunculata C. K. Allen 
L. metcalfiana C. K. Allen 
Sect. Cupuliformes H. P. Tsui 
L. megaphylla Hemsl. 
Sect. Daphnidium (Nees) Benth. 
L. chunii Merr. 
L. pulcherrima (Nees) Hook. f. 
L. villipes H. P. Tsui 
Sect. Lindera 
L. kariensis W. W. Sm. 
L. reflexa Hemsl. 


Sect. Palminerviae Meisn. 
L. obtusiloba Blume 


Sect. Polyadenia (Nees) Benth. 
L. communis Hemsl. 
Sect. Sphaerocarpae H. P. Tsui 
L. fruticosa Hemsl. (= L. neesiana 
(Wall. ex Nees) Kurz) 


Litsea Lam. 


Sect. Conodaphne (Blume) Benth. & 
Hook. f. 

a monopetala (Roxb.) Pers. 

. umbellata (Lour.) Merr. 


M — — 


n variabilis Hemsl. var. variabilis 
L. variabilis var. oblonga Lecomte 
Sect. Cylicodaphne (Nees) Hook. f. 


Voucher Source ETS ITS 
Wuyishan Exped. 1693 (KUN Guizhou, China AY817123 AY817113 
0601976) 
Li H.-W. 2 (HITBC) Yunnan, China AY934881 AY265399 
Li J. 2002032 (HITBC) Yunnan, China AY817130 AY817120 
Li H.-Q. 40091 (KUN 0100643) Guangxi, China AY817124 AY817114 
Li C.-Q. 3979 (IBK 00003410) Guangxi, China AY817122 AY817112 
Li H.-W. 1 UIT BC) Yunnan, China AY934880 AY 265398 
Yang G.-H. 55824 (KUN 0047252) Yunnan, China AY817127 AY817117 
Hainan Exped. 962 (IBK 00003425) Hannan, China AY817128 AY817118 
Xie L.-S. 613 (KUN 0047277) Guangxi, China AY817125 AY817115 
Sun B.-X. 0757 (KUN 0047286) Yunnan, China AY817126 AY817116 
Feng G.-M. 22638 (KUN 0047322) Yunnan, China AY817129 AY817119 
Wu C.-Y. et al. 75-1048 (KUN Tibet, China DQ120552 AY265397 
0049206) 
Li H.-W. 27 (HITBC) Yunnan, China 50120551 AY265396 
Li H.-W. 16 (HITBC) Yunnan, China DQ120553 AY265392 
Li J. 2002076 (HITBC) Yunnan, China DQ120541 DQ124264 
Li H. 11114 (KUN 0690203) Yunnan, China 50120542 DQI24265 
Li H.-W. 8 (HITBC) Yunnan, China 50120543 AY265408 
Li H.-W. 7 (HITBC) Yunnan, China AY934882 AY265406 
Li J. 2002134 (HITBC) Guangxi, China 50120547 001242006 
Li J. 2002202 (HITBC) Guangxi, China 50120548 00124207 
Li H. et al. 11876 (KUN 0762162) Yunnan, China 50120550 00124268 
Li H. et al. 15305 (KUN 0789698) Yunnan, China 50120539 50124203 
Nei M.-X. & Lai S. K. 3768 (KUN Jiangxi, China DQ120540 AY265407 
0100201) 
Sino-Amer. Exped. 1308 (KUN Hubei, China 0120546 AY265411 
0151469) 
Li H.-W. 4 (HITBC) Yunnan, China DQ120544 AY265409 
Li G. F. 63966 (KUN 0104915) Sichuan, China DQ120545 AY265410 
Li J. 2002108 (HITBC) Guangxi. China DQI2052; DQ120002 
Li H.-W. 24 (HITBC) Yunnan, China DQ120528 AY265404 
Li J. 2002142 (HITBC) Guangxi, China 0120529 09120003 
Li J. 2002115 (HITBC) Guangxi, China DQ120530 DQ120604 
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Table 1. Continued. 
Taxon 
L. acutivena Hayata 


~ 
* 


. dillenüfolia P. Y. Pai & 
P. H. Huang 
„n elongata (Nees) Hook. f. 
L. garrettti Gamble 
„ liyuyingi H. Liu 
» longistaminata (H. Liu) Kosterm. 
, panamanja (Nees) Hook. f. 


— —— H— — — 


a yaoshanensis Yen C. Yang & 
P. H. Huang 
Sect. Litsea 
L. glutinosa (Lour.) C. B. Rob. 
Sect, 
.. cubeba (Lour.) Pers. 
.. kingii Hook. f. 


. rubescens Lecomte 


~ 


— X — 


„ sericea (Wall. ex Nees) Hook. f. 
Neolitsea (Benth.) Merr. 


N. aurata (Hayata) Koidz. var. aurata 

N. aurata var. chekiangensis (Nakai) 
Yen C. Yang & P. H. Huang 

N. brassii C. K. Allen 


N. cambodiana Lecomte var. glabra 
C. K. Allen 

N. chrysotricha H. W. Li 

N. chuii Merr. 

N. confertifolia (Hemsl.) Merr. 

N. dealbata (R. Br.) Merr. 


N. homilantha C. K. Allen 

N. kwangsiensis H. Liu 

N. levinei Merr. 

N. lunglingensis H. W. Li 

N. ovatifolia Yen C. Yang & P. H. 
Huang var. ovatifolia 

V. ovatifolia var. puberula Yen C. 
Yang & P. H. Huang 

N. pallens (D. Don) Momiy. & H. 
Hara 

N. phanerophlebia Merr. 

N. pingbienensis Yen C. Yang & 
P. H. Huang 

N. pinninervis Yen C. Yang & P. H. 
Huang 

V. polycarpa H. Liu 


N. pulchella (Meisn.) Merr. 

N. sericea (Blume) Koidz. 

Neolitsea sp. 

N. sutchuanensis Y. C. Yang var. 
sutchuanensis 

N. sutchuanensis var. gongshanensis 
H. W. Li 

N. undulatifolia (H. Lév.) C. K. Allen 


Tomingodaphne (Blume) Hook. f. 


Li J. 2002199 (HITBC) 
Li H.-W. 19 (HITBC) 


Li J. 2002146 (HITBC) 
Xia Y.-M. s.n. (HITBC) 
Xia. Y.-M. s.n. (HITBC) 
Xia Y.-M. s.n. (HITBC) 
Li J. 2002028 (HITBC) 
Li J. 2002198 (HITBC) 


Li H.-W. 21 (HITBC) 


Li H.-W. 26 (HITBC) 

Li J. 2002170 (HITBC) 

Li J. 2002094 (HITBC) 

Li H. et al. 15299 (KUN 0789681) 


Li J. 2002181 (HITBC) 
Zhang S.-Y. 5482 (KUN 0162041) 


Gray, B. 03911 (KUN 0793628) 
Li X.-G. 202474 (IBK 00009945) 


Wu 5.-G. 7095 (KUN 0106438) 
Li J. 2002063 (HITBC) 

Xi X.-Y. 414 (PE 1272040) 
Gray, B. 03993 (KUN 0793630) 


Li J. 2002071 (HITBC) 
Wu S.-J. 3419 (IBK 00010186) 
Li H.-W. 29 (HITBC) 

Li J. 2002058 (HITBC) 

Wu S.-J. 3246 (IBK 00010360) 


Mao P.-Y. 03875 (KUN 0108307) 


Qinhai-Tibet Exped. 5972 (KUN 
0108358) 

Deng L. 7511 (KUN 0108338) 

Mao P.-Y. 04139 (KUN 0108220) 


Li J. 2002187 (HITBC) 


Zhou Z.-K. et al. EXLS-0252 (KUN 
0695075) 

Li J. 2002166 (HITBC) 

K. Midorikawa 2180 (KUN 0108215) 

Li J. 2002070 (HITBC) 

Zhao Z.-X. 64 (KUN 0108178) 


Feng G.-M. 6987 (KUN 0108134) 


Li J. 2002203 (HITBC) 


Source 


Guangxi, China 
Yunnan, China 


Guangxi, China 
China 
China 
China 
China 
Guangxi, China 


Yunnan, 
Yunnan, 
Yunnan, 
Yunnan, 


Yunnan, China 


Yunnan, China 
Guangxi, China 
Yunnan, China 
Yunnan, China 


Guangxi, China 
Zhejiang, China 


Queensland, 
Australia 
Guangdong, China 


Yunnan, China 
Yunnan, China 
Hunan, China 
Queensland, 
Australia 
Yunnan, China 
Hong Kong, China 
Yunnan, China 
Yunnan, China 
Hong Kong, China 


Yunnan, China 
Tibet, China 


Guangdong, China 
Yunnan, China 


Guangxi, China 
Yunnan. China 
Guangxi, China 
Honshu, Japan 
Yunnan, China 
Sichuan, China 


Yunnan, China 


Guangxi, China 


DQ120531 
DQ120532 


DQ120533 
0120534 
50120535 
50120530 
0120537 
50120538 


AY934883 
DQ120523 
0120524 


0120525 
50120520 


0120557 
0120556 


50120559 
0120500 
50120501 
0120562 
0120563 
50120504 
50120505 
50120560 
AY934884 
DQ120567 
DQ120568 
0120500 
0120570 


0120521 
0120572 


50120575 
5012057. 
90120575 
50120576 
0120581] 
DQ120577 
DQ120578 


0120570 


0120605 
41205105 


50120000 
0120007 
50120008 
50120000 
50120010 
DQ120611 


AY265403 
AY265402 
DQ120599 


DQ120600 
DQ120601 


50124270 
0124271 


0124272 
50124273 
50124274 
50124275 
50124276 
0124277 
50124278 
DQ121270 
AY265401 
DQ124280 
DQ124281 
DQ124282 
DQ124283 


0124284 
0124285 


0124280 
0124287 
50124288 
50124280 
DQ124294 
DQI24290 


DQ124291 


DQI2 1202 


Volume 95, Number 4 


Li et al. 585 


2008 Phylogenetic Relationships of Litsea and Laureae 
Table 1l. Continued. 
Taxon Voucher Source ETS ITS 
N. wushanica var. pubens Yen C. Yang Liu L.-H. 15149 (KUN 0162057) Hunan, China DQ120580 DQ124293 
& P. H. Huang 
Parasassafras D. G. Long 
P. confertiflora (Meisn.) D. G. Long Qian Y.-Y. 682 (KUN 0104558) Yunnan, China AY934885 AY265395 
Sinosassafras H. W. Li 
5. flavinervia (C. K. Allen) H. W. Li Liu Y.-H. s.n. (HITBC) Yunnan, China AY934886 AY940451 
Outgroup taxa 
Cinnamomum pittosporoides Hand.-Mazz. Li H. 5252 (KUN 0108156) Yunnan, China DQ120554 DQ124269 
Sassafras tzumu (Hemsl.) Hemsl. Li H.-W. 15 (HITBC) Yunnan, China DQ120555 AY265391 
Umbellularia californica (Hook. & van der Werff s.n. (MO) North America DQ120556 AY265393 


Arn.) Nutt. 


implemented in PAUP*, and because the PHT was 
non-significant, the data sets were combined for all 
subsequent analyses. Most parsimonious trees were 
obtained from 10,000 replicates of random taxon 
addition using equally weighted (EW) characters 
(Fitch, 1971) and tree bisection-reconnection (TBR) 
(MULPARS off, 
swapping on the shortest trees from this analysis with 
MULPARS on. These trees were then used to re- 


weight the characters according to the best fit of their 


branch swapping followed by 


rescaled consistency indices (Farris, 1989). New 
searches were performed with 1000 replicates using 
successive weighting (SW) until equilibrium was 
reached (Farris, 1969). Clade support was estimated 
using bootstrap resampling (Felsenstein, 1985), with 
1000 replicates, TBR, EW, and MP optimality criteria 
performed on the combined weighted ITS and ETS 


Table 2. 


generic and sectional level in Laureae. 


Habit: evergreen (0), deciduous (1) 


data with resampling using all characters equally, 
regardless of weight. 

The 19 morphological characters were analyzed 
using unweighted parsimony, with all characters 
treated as unordered. Tree search was performed with 
1000 replicates of random taxon addition and TBR 
branch swapping (MULPARS on) in PAUP* version 
LOLbLO (Swofford, 1999). The morphological data 
were then combined with the ITS and ETS data and 
analyzed using the same settings that were used in the 
morphological data set. Clade support was estimated 
with bootstrap resampling (Felsenstein, 1985) for 
1000 replicates with TBR, EW, and MP optimality. 
Morphological character state changes were then 
plotted on one of the resulting most parsimonious 
combined analysis trees using MacClade (Maddison & 


Maddison, 2000). 


Morphological characters and character states traditionally considered to be taxonomically important at the 


l 
2. Leaf arrangement: alternate along stems (0), alternate and crowded at branchlet apices (I), verticillate or subverticillate (2) 


3. Leaf venation: pinninerved (0), triplinerved (1), trinerved (2) 


4. Inflorescence type: thyrsoid, without vegetative terminal bud in the main axis (0), short shoot (brachyblast) with 


vegetative terminal bud in the main axis (1) 


5. Inflorescence: terminal or subterminal (0), axillary (1) 


6. Inflorescence arrangement: panicle (0), raceme (1), fasciculate clustered (2) 


7. Inflorescence: sessile (0), stipitate (1) 

8. Flower number per inflorescence: >1(0), 1(1) 
9. Involueres: absent (0), present (1) 
10. Involucres: large (0), minute (1) 
ll. Involucres: imbricate (0), decussate (1) 
12. Involucres: early deciduous (0), persistent (I) 
13. Flower sex: bisexual (0), unisexual (1) 

14. Basie floral number: dimerous (0), trimerous (1) 


16. Anthers: two-locular (0), four-locular (1) 


Perianth segment: present, perfect (0), imperfect, absent or early deciduous (1) 


17. Pollen sacs of the third whorl: latrorse (0), introrse (I), extrorse (2) 


18. Fruit shape: globular or oblate (0), ovoid or ellipsoid (1) 


19. Fruit cupule shape: flat or discoid (0), cup-shaped (1) 
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Bayesian phylogeny reconstruction of the combined 
data was performed with MrBayes 3.0b4 (Huelsenbeck 
& Ronquist, 2001: Huelsenbeck et al., 2002). The 
program Modeltest (Posada & Crandall, 1998) was used 
to test which of the 56 predicted. models of DNA 
substitution best fit the data. Modeltest indicated that 
the FrN + 1 + F model with rate heterogeneity and 
among-sile rate variation from Tamura and Nei (1993) 
best fit the two nrDNA regions, and the ML parameters 
in MrBaves (MB) were "Ist nst = 67 


inveamma." The Markov chain Monte Carlo process was 


and "rales = 


set so that four chains ran simultaneously for 500.000 
generations, with trees sampled every 100 generations. 
giving a total of 5000 trees in the initial sample. 
Likelihood value plots for the four chains showed that 
stationarity had occurred by the 600th tree. Therefore. 
the first 600 trees were discarded as “burn in.” and the 
posterior probabilities of the phylogeny and its branches 
were determined from the remaining 4400 trees. 

ME analysis was also performed using neighbor 
joining on the unweighted data set using MEGA 3.1 
(Kumar et al., 2004). with bootstrap support calculat- 
ed on 10.000 replicates. 

ML analysis of the combined, unweighted data was 
performed using the DNAML option in DAMBE 
version 4.13 (Xia, 2000; Xia & Xie. 2001). 


RESULTS 
SEQUENCE CHARACTERISTICS 


[TS regions are quite variable between distantly 
related taxa. so only regions that could be aligned 
unequivocally were used in this analysis, making our 
phylogenetic estimates conservative. For the ingroup 
laxa, the length of the ITS regions. including the 5.95 
region, ranged from 568 to 627 bp, and the alignment 
generated a data set of 689 characters of which 156 
(22.04%) were parsimony informative. G plus C 
content ranged from 64.61% to 74.19%. 

Compared to ITS, the approximately 400 bp of the 
3' end of the ETS between the 185 subunit and the 
ETS-1 internal primer were relatively easy to amplify 
and sequence. For the ingroup taxa. the ETS region 
varied in length from 350 to 393 bp. and alignment 
resulted in a matrix of 393 characters of which 99 
(25.19%) were parsimony informative, with G plus C 
46.56% to 55.23%. All 


sequences were submitted to GenBank (Table 1). 


content ranging from 


and the data matrices for both sequenced regions are 
available from the primary author upon request. 


PHYLOGENETIC ANALYSIS 


SW parsimony. Results of phylogenetic analyses 


of the ITS and ETS regions performed separately 


showed no hard incongruences (i.e.. there were no 
o 
contradictory clades supported by bootstrap greater 
than 60%: data not shown). The PIHT (Farris et al.. 
1995) indicated that the two regions were congruent 
D D 
(P^ « (LI: see 1996; 
Cunningham, 1997: Farris et al., 2000). Because of 


discussions. in Sullivan. 
this and. because the ETS and ITS regions. occur 


within the same transcriptional unit and show 
evidence of a similar and interdependent role in the 
maturation of ribosomal RNAs (Good et al.. 1997), we 
will only present results for the combined analyses. 
The EW analysis produced 30.600 trees (length [L] 
= 1171. consistency index [CI] — 0.306. retention 


index [RI]. = 


virtually all minor branches collapsed under strict 


0.589). but because all main and 


consensus. SW was used to try to stabilize the tree 
topology. SW reduced this to 12 trees (L = 211.5, CI 
= 0.576. RI = 
resolved strict consensus tree (Fig. |). 

The SW analysis of the nrDNA data produced a 


series of clades within a monophyletic core Laureae 


0.834). producing an almost fully 


(93% bootstrap support). These clades are referred to 
informally as the Neolitsea—Actinodaphne, and Litsea, 
Lindera, and Aperula clades. These sat above a basal 
grade of Lindera obtusiloba Blume (sect. Palminer- 
piae) and then Lindera communis Hemsl., (sect. 
Polyadenia). The Aperula clade contained the three 
sampled species of Lindera sect. Aperula (I. latifolia, 
L. longipedunculata C. K. Allen. and £L. metcalfiana €. 
K. Allen) grouped in a terminal pairing with Litsea 
cubeba and L. kingii Hook. f. (both section Tomingo- 
daphne) above a subclade of Lindera. megaphylla 
(ect. Cupuliformes) and Actinodaphne forrestii (C. K. 
Allen) Kosterm. Although the Aperula clade as a 
whole was unsupported (< 50% bootstrap), all the 
branches within it had moderate (> 70%) to strong 
support (> 90%). 

Siling above Sinosassafras and Parasassafras was 
the Lindera clade (75% bootstrap support). with two 
subelades. The first of these represented Lindera sect. 
Lindera (2 spp.). plus Litsea species from sections 
Tomingodaphne and Conodaphne (part). Sister to this 
was a subclade consisting of L. fruticosa Hemsl. (sect. 
Sphaerocarpae). leadaphne caudata (Nees) H. W. Li. 
and the three species of Lindera sect. Daphnidium. 


again wilh most terminal branches moderately 
supported, 


The Litsea clade consisted of a terminal lineage 
representing members of sections Litsea. Conodaphne 
(part). and Cylreodaphne, but also including Dodeca- 
denia grandiflora Nees, all sitting above Laurus 
Allen, 


The clade included within it three separate lineages. 


nobilis U, and Actinodaphne lecomtei C. K. 


Section Litsea and Litsea monopetala (Roxb.) Pers. 


(sect. Conodaphne) formed a strongly supported pair 
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(99%), sister to a subclade to three species of Litsea 
sect. Cylicodaphne: L. dillenifolia P. Y. Pai & P. H. 
Huang, I. garrettii Gamble, and L. panamanja (Nees) 
Hook. f. (here called Cylicodaphne V), and then to an 
L. variabilis Hemsl. var. variabilis and var. oblonga 
Lecomte (sect. Conodaphne) subclade, all with 
moderate to strong bootstrap support. The other 
branch in the Litsea clade was the Cylicodaphne H 
subclade, representing the remainder of Litsea sect. 
Cylicodaphne but including an embedded Dodecade- 
nia grandiflora, agam with most branches showing at 
least moderate support. 

Within the Neolitsea—Actinodaphne clade, Actino- 
daphne (except A. forrestii and A. lecomtei) formed a 
basal grade to a well-supported Neolitsea (89%). with 
the latter divided into two subclades (Neolitsea | and 
Neolitsea M) above N. chrysotricha H. W. Li and N. 
pallens (D. Don) Momiy. & H. Hara. Neolitsea | 
showed little clear support for the internal branches. 
Neolitsea |l showed 


whereas all the branches in 


bootstrap support > 50%. 


Bayesian analysis. Bayesian analysis of the 


unweighted Laureae [TS + ETS showed moderate 
relationship resolution (Fig. 2), with the Bayesian tree 
corresponding well with much of the SW tree in terms 
of recovered major lineages. Although not as well 
resolved as the SW topology, terminal SW clades with 
high bootstrap support were also present in the 
Bayesian topology with strong posterior probability 
support, and both analyses included the Neolitsea— 
Actinodaphne clade and many of the major SW 
subclades. 

Nevertheless. there were differences between the 
results for the two approaches. In the Bayesian tree, 
the Lindera clade was not recovered. with section 
Daphnidium falling instead as part of a polytomy 
separate from the remainder. Similarly, Litsea gluti- 
nosa (Lour.) C. B. Rob. (sect. Litsea and the type 
species for the genus) and L. monopetala (sect. 
Conodaphne) were separated from the rest of the 
Litsea clade seen in the SW analysis. 


ME and ML analyses. 
tree (Fig. 3) was well resolved and similar in major 


The ME (neighbor joining) 


clade structure to the SW and Bayesian cladograms. 
Neolitsea 1 and II (albeit the latter reduced) were 
again terminal above an Actinodaphne grade, although 
here A. paotingensis Yen C. Yang & P. H. Huang was 
embedded in a basal Neolitsea subclade with N. 
brassii C. K. Allen, N. wushanica var. pubens Yen C. 
Yang & P. H. Huang, and N. sutchuanensis Y. C. 
Yang. Laurus was basal to the Actinodaphne—Neolitsea 
clade. and sister to this group was a clade consisting 
of part of the MP and MB Aperula clade and an 
expanded Litsea clade including the Litsea and 


subelades. as well as 


and II 


Dodecadenia and lteadaphne. Below this was a 


Cylicodaphne | 


lineage containing the Lindera, L. fruticosa, and 


Daphnidium  subclades, and then a polytomy 
consisting of L. communis paired with Sinosassafras, 
and an 


"Aperula I” clade of Tomingodaphne, 


Cupiliformes, and Actinodaphne forrestii. Lindera 
obtusiloba and then Parasassafras were basal to the 
remainder of Laureae. 

Unlike most of the other analyses, in the ML tree 
(Fig. J). Lindera obtusiloba and L. communis were 
paired and sister to the Aperula clade, with this 
lineage placed above Laurus and below the Actino- 
daphne—Neolitsea clade. Neolitsea J and II and the 
Actinodaphne grade were all present, although A. 
paotingensis was placed inside Neolitsea I, similar to 
the ME tree, and a subclade of N. aurata. (Hayata) 
Koidz. and N. kwangsiensis H. Liu was basal to 
Veolitsea as a whole. Below all of these was a major 
lineage consisting of a Lindera clade, L. fruticosa, and 
the Daphnidium clade (with /teadaphne basal) and a 
Litsea clade (including Litsea and the Cylicodaphne | 
and H subclades). Actinodaphne lecomtei was placed 


well inside the Actinodaphne grade, whereas A. 


forrestii was still sister to Lindera megaphylla, but 


this latter pair was placed between Sinosassafras and 


Parasassafras at the base of the Laureae. 


Combined molecular and morphological analyses. 
Analysis of the morphological matrix by itself resulted 
in 317 equally parsimonious trees (L = 72, CI = 
0.3056, RI = 0.8214), but these collapsed completely 
to an unresolved polytomy under strict consensus. 
When the morphological and molecular data were 
combined, two trees of 1271 steps (CI — 0.3021, RI — 
0.6086) resulted, one of which is shown in Figure 5. 
This tree also recovered the Neolitsea—Actinodaphne 
clade, Lindera clade, and Litsea clade, which were 
seen in the molecular analysis, but the Aperula clade 
was now split into two separate entities, with Aperula | 
representing Lindera sect. Aperula s. str. and Aperula 
Il representing Litsea cubeba and L. kingii (sect. 


Tomingodaphne), Lindera megaphylla (sect. Cupuli- 


formes), and Actinodaphne forrestii. The inclusion of 


the morphological data changed the bootstrap support 
for most clades, generally lowering it from the SW 
analysis, although there was still reasonable bootstrap 
support for the terminal branches and still little or no 
deep-branch support. 

When the morphological character state changes 
are plotted on the combined analysis tree (Fig. 5), the 
unique synapomorphy for the Laureae is the presence 
of inflorescences with short shoots (brachyblasts) with 
a vegetative terminal bud on the main axis. Within the 
Laureae, the Actinodaphne—Neolitsea—Parasassafras— 
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Table 3. Data matrix of important morphological characters in Laureae for species used in the molecular analyses. 


— 
© 
KA 
— 
= 
NO 
— 
LA 
p 
VS 


15 16 17 18 19 


Taxa 


Actinodaphne cupularis 
Actinodaphne forrestii 
Actinodaphne henryi 
Actinodaphne kweichowensis 
Actinodaphne lecomtei 
Actinodaphne obovata 
Actinodaphne omeiensis 


Actinodaphne paotingensis 
Actinodaphne pilosa 


Actinodaphne trichocarpa 


eS GOGO OOOO O GO 
SOS SS 


Actinodaphne tsati 
Cinnamomum pittosporoides 
Dodecadenia grandiflora 
Iteadaphne caudata 


Laurus nobilis 


C 


Lindera chunti 
Lindera communis 


Lindera fruticosa 


à 


Linder kariensis 

Lindera latifolia 

Lindera longipedunculata 
Lindera megaphylla 
Lindera metcalftana 
Lindera obtusiloba 
Lindera pulcherrima 
Lindera reflexa 


GEN Gel COME CO (o 


Lindera villipes 


Litsea aculivena 
Litsea cubeba 
Litsea dillentifolia 
Litsea elongata 
Litsea garretti 
Litsea glulinosa 


e OOGO a H E A E ta kär 


Litsea kingii 
Litsea liyuyingi 
Litsea longistaminata 


Litsea monopetala 


„ 


Litsea panamanja 


Litsea rubescens 


el 


Litsea sericea 
Litsea umbellata 
Litsea variabilis var. variabilis 


Litsea variabilis var. oblonga 


Litsea yaoshanensis 

Neolitsea aurata var. aurata 
Neolitsea aurata var. chekiangensis 
Neolitsea brassii 

Neolitsea cambodiana var. glabra 


Neolitsea chrysotricha 


Neolitsea chuii 
Neolitsea confertifolia 
Neolitsea dealbata 


EEE 


Neolusea homilantha 
Neolitsea kwangsiensis 


Neolitsea levinei 
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Neolitsea lunglingensis 
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Continued. 


Table 3. 


Taxa 


Neolitsea ovatifolia var. ovatifolia 
Neolitsea ovatifolia var. puberula 
Neolitsea pallens 

Neolitsea phanerophlebia 
Neolitsea pingbienensis 

Neolitsea pinninervis 

Neolitsea polycarpa 

Neolitsea pulchella 

Neolitsea sericea 

Neolitsea sp. 


OO C) C9 OO E e en E em De Ween D 
O0 0 
E= EEE O EEE 
GOS OO) 'O Q(VO OO OO. 
E= EEE 
O O D D DO 


Neolitsea sutchuanensis var. 
sutchuanensis 


CH 
CH 
— 
CH 
— 
No 


Neolitsea sutchuanensis var. 
gongshanensis 

Neolttsea undulatifolia 

Neolitsea wushanica var. pubens 

Parasassafras confertiflora 

Sassafras tzumu 

Sinosassafras flavinervia 


Shem Ren eme 
QUO C OR 
0828 
S0 
20 
E 


Umbellularia californica 


Sinosassafras clade is united by the presence of 
imbricate, early deciduous involucral bracts, although 
Neolitsea itself has decussate persistent bracts as one 
of its synapomorphies. Parasassafras and Sinosassa- 
fras are supported as a pair by the possession of 
triplinerved leaves, minute involucral bracts, and 
latrorse pollen sacs for the third whorl, whereas the 
monophyly of the Neolitsea—Actinodaphne clade is 
supported by the synapomorphies of verticillate or 
subverticillate leaves and thyrsoid inflorescences that 
lack a terminal bud. Although there are few characters 
supporting the main branches within the Actinodaphne 
grade, Neolitsea is well supported, with triplinerved, 
alternate, crowded leaves at branch apices, decussate 
Within 


Neolitsea, there were no morphological synapomor- 


persistent bracts, and dimerous flowers. 


phies supporting clades I and II, but within I, the 
species above N. ovatifolia Yen C. Yang & P. H. 
Huang var. ovatifolia (Fig. 5la) have verticillate or 
subverticillate leaves as a synapomorphy (albeit with 
some reversals), whereas in clade II there are two 
lineages: Ila, with ovoid to ellipsoid fruits, and IIb, 
defined by alternate leaves (again with a reversal in N. 
brassit). 

In contrast, there were relatively no morphological 
synapomorphies supporting the Lindera or Aperula | or 
II clades, although the Laurus—Litsea clade had 
racemose inflorescences, and the Litsea clade was 
defined by four-locular anthers and the L. glutinosa— 
L. monopetala pair by the possession of imperfect, 
absent or early deciduous perianth segments. Within 
the Litsea clade, the Cylicodaphne II subclade was 
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characterized by cup-shaped cupules and united with 
Dodecadenia on fasciculate, clustered inflorescences 
and ovoid/ellipsoid fruits. Similarly, within the 
Lindera clade, section Daphnidium was united by 
trinerved leaves and ovoid-ellipsoid fruits. 


DISCUSSION 
PHYLOGENETIC UTILITY OF ETS IN THE LAUREAE 


Several studies have shown that greater resolution 
and support for phylogenetic estimation are achieved 
by increasing character number and/or taxon repre- 
sentation (Graybeal, 1998; Hillis, 1998; Soltis et al., 
1998; Bremer et al., 1999), and our study supports the 
importance of adding data from both more taxa and 
more sequence regions to help resolve issues in 
Laureae. The usefulness of the ETS region in 
molecular systematics has been suggested previously 
(Baldwin & Markos, 1998; Bena et al., 1998; 
Clevinger & Panero, 2000; Linder et al., 2000), and 
although the sequenced ETS segment is much shorter 
than the segment in the ITS (393 bp vs. 689 aligned 
bp) in Laureae, it nevertheless produces a slightly 
higher percentage of informative sites (25.19%) than 
that in the ITS (22.64%). The PHT for the two data 
sets showed congruence, and the combined analysis of 
the ITS and ETS sequences provided greater resolu- 
tion and increased support for the relationships than 
either sequence by itself (trees not shown). This is 
consistent with the results of other combined ITS/ETS 
investigations (e.g., Li et al., 2002; Becerra, 2003; Lee 
et al., 2003; Morgan, 2003; Saar et al., 2003; Urbatsch 
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Neolitsea/Actinodaphne Clade 
52 


IN 


Litsea Clade p 99 


Lindera Clade "2 9 


Aperula Clade op 


sequence data. Bootstrap values greater than 50% are indicated on branches. 


et ul., 2003: Plovanich & Panero. 2004: Roalson & 


Friar, 2004). 


RELATIONSHIPS AMONG MAJOR CLADES 


Laureae are a well-supported monophyletic group 


based on our molecular data. and they share a 


Neolitsea cambodiana var glabra 
Neolitsea ovatifolia var puberula 
Neoltsea dealbata 

Neolitsea phanerophiebia 
VDH pukhella 

Neoltsea sericea 

Neolitsea levine; 

Neoltsea ovatifolia var ovatifolia 
Neoltsea brassii 

Neolitsea aurata var. aurata 


Neoltsea aurata var. chekiangensis 


Neolitsea kwangsiensis 
Neoltsea funglingensis 
Neoltsea confertifolia 
Neolttsea undutatfolia 
Neolitsea chuii 


Neolitsea sutchuanensis var sutchuanensis 
Neolftsea sutchuanensis var gongshanensis 


Neolitsea wushanica var pubens 
Neolitsea pingbienensis 
Neolitsea polycarpa 
Neolitsea homilantha 
Neolitsea sp 

Neolftsea pinninervis 
Neolitsea chrysotricha 
Neoltsea pallens 
Actinodaphne paotingens!s 
Actinodaphne cupularis 
Actinodaphne kweichowensis 
Actinodaphne omelensis 
Actinodaphne trichocarpa 
Actinodaphne pilosa 
Actinodaphne tsaii 
Actinodaphne henryi 
Actinodaphne obovata 
Litsea garrettii 

Litsea panamanja 

Litsea dillenifolia 

Litsea variabilis var. variabilis 
Litsea variabilis var. oblonga 
Litsea glutinosa 

Litsea monopetala 

Litsea acutivena 

Litsea , aoshanensis 

Litsea elongata 
Dodecadenia grandifiora 
Litsea liyuyingi 

Litsea longistaminata 
Actinodaphne lecomter 
Laurus nobilis 

Lindera pulcherrima 
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morphological synapomorphy of pseudo-umbellate, 


rarely racemose or paniculate involucrate inflores- 


cences. The results of our analyses are. in general, 


similar to the ITS + matK phylogeny of Li et al. 


(2004). and the four clades of Laureae and other 


subgroups recognized in that study more or less 


correspond to clades recovered here. Increased taxon 
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Phylogenetic majority rule tree for Laureae from Bayesian inference analysis of the unweighted combined ITS 


and ETS sequence data showing posterior probabilities on branches. 


sampling as well as the additional sequence data from 
the ETS region and morphology have largely improved 
bootstrap, but mainly in more derived regions, with 
This 


suggests that although main lineages are being 


deep branch relationships still unresolved. 


identified and supported, the positions of many of 


these clades are unstable under different analysis 


procedures (e.g., Lindera obtusiloba basal to the 
Laureae; the isolated, generally basal position of L. 
communis; and the clades Lindera and Aperula less 
derived than the Litsea, Laurus, and Actinodaphne— 
Neolitsea clades). The precise relationships between 
major clades are, thus, still uncertain, suggesting that 


the sequence regions are evolving too rapidly to provide 
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Figure 3. 


Minimum evolution (neighbor joining) tree for Laureae of the unweighted combined [TS and ETS sequence 


data. Bootstrap values (10.000 replicates) receiving support greater than 50% are indicated on branches. 


deep-branch resolution, and sequences from slower- 
evolving gene regions may be needed to stabilize and 
support relationships between the major clades. 


However. despite this. caveat, our study includes 


representative taxa from all genera and most sections of 


Laureae, and the resulting clades support the hypothesis 


of Li et al. (2004) that inflorescence features and 
onlogeny are important for helping to understand 
evolution and improve classification within the Laureae. 

The moderately supported association. between 
Sinosassafras flavinerria and Lindera communis (sect. 


Polyadenia) in the MB and ME analyses also suggests 
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Figure 4. 


that the former may not be as closely related. to 


Parasassafras as previously suggested (Rohwer, 


1993), 


successively as basal to the Lindera clade in the SW 


although these two genera were placed 
analysis (albeit without bootstrap support) and formed 


a clade in the combined molecular/morphology 
analysis, with three morphological synapomorphies. 
Sinosassafras flavinerria has a single or two(to three) 


pseudo-umbels clustered in the leaf axils and always 
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bears several pseudo-umbels on axillary short shoots, 
whereas L. communis only has a single or rarely two(to 
three) pseudo-umbels clustered in leaf axils. Never- 
theless, in both taxa, the terminal buds on the axillary 
short shoots are poorly developed or reduced. 

The Aperula M clade contained two subgroups (none 
of which actually belong to Lindera sect. Aperula s. 
str.): (1) Litsea cubeba and L. kingii; and (2) Lindera 
megaphylla and Actinodaphne forrestii, and three of 
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Figure 5. One of the two trees resulting for Laureae from unweighted combined ITS, ETS. and morphological data sets. Morphological character state changes are also plotted. Numbers 
above branches are characters: those below are character states from Table 2. Filled circles are unique state changes: open circles represent homoplasious changes. Bootstrap values greater than 


50% are indicated in brackets below branches. 
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these (Litsea cubeba, Lindera megaphylla, and 4. 
forrestii) also formed a clade in the study by Li et al. 
(2004). Litsea cubeba and L. kingit (sect. Tomingo- 
daphne) are united by being deciduous, with naked 
terminal buds, and differ mainly in the absence of leaf 
pubescence in £L. kingii (Long, 1984). Li et al. (2004) 
also noted that there was a close relationship between 
Lindera megaphylla (sect. Cupuliformes) and A. 
forrestii, and that there were micromorphological 
grounds for the splitting of Actinodaphne s.l. This is 
further supported by a preliminary molecular study of 
Actinodaphne, which found the genus to be polyphy- 
letic within Laureae (Li et al., 2006). Although Li 
(1985) and Tsui (1987) had previously regarded 
Lindera sect. Cupuliformes as being possibly related 
to Litsea sect. Cylicodaphne because of their similar 
fruit cupules (despite belonging to different genera), 
this was not supported by our study, in which section 
Cylicodaphne itself was polyphyletic. Furthermore. 
Lindera megaphylla has large leaves aggregated near 


the top of branchleis (an unusual character in 
Lindera), making it morphologically more similar to 
some Actinodaphne species. Lindera megaphylla also 
bears a pair of pseudo-umbels on each side of an 
axillary short shoot with a vegetative terminal bud, 
whereas A, forrestii has several sessile pseudo-umbels 
clustered on an axillary short shoot that also produces 
a vegetative terminal bud. This suggests that L. 
megaphylla and A. forrestii may share a common 
inflorescence ontogeny and could help to explain why 
A. forrestii was separated from other Actinodaphne 
species in the analyses of Li et al. (2004, 2006) and 
the current study. 

The largely basal and isolated position of Lindera 


obtusiloba (sect. Palminerviae) is unusual. This 
species is deciduous, and its pseudo-umbels are 
borne in an axillary, mixed bud (leaves and 


inflorescences together) covered by scales. This kind 
of mixed bud also occurs in Sassafras. According to 
this character and the later appearance of L. 
obtusiloba in the fossil record, Tsui (1987) suggested 
that Lindera evolved from Sassafras during the 
Miocene, possibly in response to the onset of cooler, 
more seasonal climates, although this assumes the 
correct assignment of fossils to extant genera, and 
generic definition based on vegetative anatomical and 
morphological features, which is an area of ongoing 
research. 

Neolitsea consistently formed a terminal monophy- 
letic lineage within the Neolitsea—Actinodaphne clade, 
agreeing with the matK analysis of Li et al. (2004). 
The present study (with many more taxa) also shows 
that Neolitsea as a clade is defined by dimerous 
flowers (apparently a reduction from a trimerous 


ancestral condition). clustered/verticillate leaf ar- 
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rangement, triplinerved venation, and decussate, 
persistent involucral bracts, and that splitting the 
genus just on leaf venation differences, as suggested 
as a possibility by Li et al. (2004), is not warranted. 

Although two main subelades (Neolitsea I and If) 
were found in all the analyses, the precise composi- 
lion of these varied, with some taxa moving between 
them or to a basal grade position depending on the 
analysis used, possibly due to what is still a relatively 
small sample size for such a large genus. Similarly, 
the absence of previously defined morphology-based 
sections within Neolitsea makes “representative” 
sampling more difficult. As they currently stand, 
neither Neolitsea | nor H has definitive synapomor- 
phies, and they do not differ consistently in leaf 
venation or inflorescence features. Leaf arrangement 
seems to be important for two of the larger subclades, 
with the species in Ia being verticillate, and those of 
IIb being alternate. Similarly, Ha above N. penninervis 
shows ovoid-ellipsoid fruit, although this feature also 
occurs in some Neolitsea species from subclades I and 
IIb. Accordingly, the composition of the clades within 
Veolitsea and the nature of any supporting morpho- 
logical characteristics are the focus of ongoing 
research. 

Similarly, the inclusion of Actinodaphne paotingen- 
sis inside Neolitsea by the ME and ML analyses 
warrants study. This placement may be related to the 
instability of the number of floral parts, as there can 
be six to eight perianth lobes and nine to 15 fertile 
stamens in this species. 

The combined molecular and morphological phy- 
logeny indicates that Neolitsea is terminal above 
Actinodaphne, despite previous morphological studies 
suggesting that Neolitsea is closest to Litsea (Koster- 
mans, 1957; Hyland, 1989; Rohwer, 1993). Van der 
Werff (2001) that the 
Actinodaphne is unlike that of other Laureae, and 


noticed inflorescence of 
several species (e.g.. A. pilosa (Lour.) Merr., A. henryi 


Gamble. and A. obovata (Nees) Blume) have a 
paniculate or racemose inflorescence enclosed by 
imbricate, early deciduous involucral bracts. In our 
study, the Neolitsea—Actinodaphne clade was defined 
by the possession of verticillate leaves (later becoming 
terminally clustered in Neolitsea) and thyrsoid 
inflorescences lacking a terminal bud. This latter 
different 


brachyblast-type short shoot seen in Laurus or the 


inflorescence feature is very from the 
Litsea, Lindera, and Aperula clades, which display 
instead a pseudo-umbellate inflorescence with a 
bud 
deciduous involucral bracts that enclose each pseu- 
do-umbel. Rohwer (1993: fig. 87C-E) suggested that 
reduction from paniculate or racemose inflorescences 


vegetative terminal and decussate, tardily 


led to the sessile pseudo-umbellate inflorescences 
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seen in Veolitsea, and this agrees with its position as a 
derived terminal above Actinodaphne. Furthermore. 
the clustered pseudo-umbels seen in Actinodaphne 
and Neolitsea differ from those of Litsea or Lindera. as 
the former are thyrsoid, lack terminal buds, and are 
clustered in leaf axils, whereas the latter are arranged 
along leafless short shoots and bear a vegetative 
terminal bud. 

Within Neolitsea, the separation of N. ovatifolia var. 
ovatifolia from variety puberula Yen C. Yang & P. H. 
Huang was unexpected. It may just reflect the 
widespread distribution and/or regional differentiation 
within this species, but certainly warrants further study 
to clarify species and varietal limits in this taxon. 

Li (1985) and Tsui (1987) suggested an evolution- 
ary series for the inflorescences in Litsea and Lindera. 
and our results concur with their hypotheses. The 
flowers occur in pseudo-umbels enclosed by decus- 
sate involucral bracts and are arranged along a 


leafless axillary short shoot with a terminal bud that 


can grow into a normal vegetative branch after 


flowering. A raceme bearing pseudo-umbels arising 


from normal growth of the peduncles and growth of the 


internodes in the short shoot was a synapomorphy for 


the Laurus—Litsea—Lindera—Aperula | clade, with 


reduction in several distal branches to create 


fasciculate clusters. In other basal groups. the 


peduncles and internodes of the short shoot are 
frequently reduced, so that the flowers are arranged in 
spikelike pseudo-umbels (the sessile inflorescence 
plesiomorphic in the 


condition appears to be 


Laureae). In /teadaphne and Dodecadenia. although 


they follow this basic racemose pattern, the number of 


flowers per involucre or pseudo-umbel is reduced to 


one. In both cases, they seem to represent reduced 


members of otherwise pseudo-umbellate clades, 
suggesting that the pattern is convergent, 
Previous morphology-based studies recognized 


considerable variability in both Litsea and Lindera 
and variously subdivided them into sections (Ben- 
tham, 1880: Hooker. 1890: Li et al., 1984; Tsui. 
1987). Hooker (1890) and Li et al. (1984) recognized 
four 


sections within Litsea (Litsea, Conodaphne, 


Cylicodaphne, and Tomingodaphne) based on habit. 


floral 


cupules. Lindera was similarly divided into eight 


leaves. characters. inflorescences. and fruit 
sections (Lindera, Sphaerocarpae. Palminerviae, Aper- 
ula, Cupuliformes, Daphnidium, Polyadenia, and 
Uniumbellatae) by Tsui (1987). Although traditional 
generic delimitations based on two-celled versus four- 
celled anthers were not supported by our study, 
several monophyletic subclades are evident that do 
correspond, in part, lo some of these previously 
recognized sections, and help to shed light on their 


phylogenetic relationships, as well as character 


evolution. within the clades. For example, the Litsea 
clade was synapomorphic for four-locular anthers, 
albeit with /teadaphne embedded within it as a 
reversal. 

Within the Litsea clade, L. glutinosa (sect. Litsea 
and type species for the genus) and 4. monopetala 
(sect. Conodaphne) were consistently resolved as 
sister taxa. with moderate support in the ME and 
Bayesian analyses. Although clustered on both 
molecular as well as combined data, the pair is 
characterized by a lack or incompleteness of the 
perianth lobes (absent in I. glutinosa and early 
deciduous in £L. monopetala) (How, 1956). 

The Cylicodaphne I and H subelades, correspond- 
ing to members of Litsea sect. Cyclicodaphne s.l., were 
well supported as separate in most of the analyses. 
and. even in the Bayesian analysis where they formed 
a single clade, Cylicodaphne | was supported and 
terminal above an unsupported Cylicodaphne H grade. 
Although I. variabilis was placed in section. Con- 
(1984). it was treated as a 


member of section Cylicodaphne by Hooker (1890) 


odaphne by Li et al. 


and our results support its return to that section. in 
Cylicodaphne I. Within Cylicodaphne J, although some 
species have eight perianth lobes. they all have 
pseudo-umbellate racemose inflorescences arranged 
along leafless short shoots, a feature that they share 
and 


subelade (at least in 


with the Litsea part) 
lteadaphne, although there were no morphological 
synapomorphies to define Cylicodaphne L. In Cylico- 
daphne II. the species have one to several clustered 
pseudo-umbels on an axillary short shoot and ovoid- 
elliptical fruits, supporting the inclusion of Dodeca- 
denia. 

Lindera fruticosa (sect. Sphaerocarpae) was related 
lo Headaphne and section. Daphnidium in the SW 
analysis, despite its possession of an umbel with an 
elongate peduncle and its deciduous habit. Neverthe- 
less, it was placed as a less-derived member of the 
Lindera clade in the ME and ML analyses. with which 
it shares a deciduous habit and possession of a single 
or few pseudo-umbels clustered on the short shoot. Its 
hasal position may be related to the intermediate 
situation of leaf venation, as it can possess. both 
triplinerved and pinninerved leaves. However, as 
none of our analyses showed bootstrap support for the 
position of this taxon, its associations must be 
regarded with caution for the present. pending further 
studies. 

The Lindera clade consists of species from section 
Lindera and members of Litsea sections Conodaphne 
and Tomingodaphne. In deciduous habit, Lindera 
kariensis W. W. Sm. (sect. Lindera) is very like Litsea 
sericea (Wall. ex Nees) Hook. f. (sect. Tomingo- 
daphne), except for the difference in number of anther 
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cells. Nevertheless, the placement of the evergreen 
Litsea umbellata (Lour.) Merr. (sect. Conodaphne) 
within the group is unusual, as all other members of 
this clade are deciduous. Litsea umbellata also 
appears unusual within the clade, as its cupule bears 
persistent tepals. Nevertheless, it is possible thal 
deciduous species that bear cupules without persis- 
tent tepals may have evolved from evergreen ancestors 
with persistent cupule tepals. and in section Con- 
odaphne s. str., there are species with and without 
persistent cupule tepals. Similarly, section Tomingo- 
daphne contains two vegetative terminal shoot bud 
forms in the axillary short shoots: naked versus scaly 
(Hooker, 1890). 

Rohwer (1993) and van der Werff and Richter 
(1996) recognized a large portion of what we now 
regard as the Laureae, including Actinodaphne, Do- 
decadenia, lteadaphne, Lindera, Litsea, Neolitsea, 
Parasassafras, and Sinosassafras. Although Sassafras 
and Umbellularia were also included, they should be 
removed based on the results of Chanderbali et al. 
(2001), Li et al. (2004), and the present study. This is 
further supported by the possession of racemose 
inflorescences in Sassafras and the presence of 
bisexual flowers with extrorse anther cells in the third 
whorl in Umbellularia. 

In our study, we sampled species representing 
all genera and most sections in the core Laureae, 
and it is clear thal a major disparity exists between 
our molecular phylogenetic results and more tradi- 
tional morphology-based taxonomic concepts of ge- 
neric and infrageneric classification in the tribe. For 
example, both current and earlier analyses (Li & 
Christophel, 2000; Li et al., 2004) confirm that using 
two- versus four-celled anthers to separate Litsea s.l. 
from Lindera s.l. results in polyphyletie or paraphy- 
letic genera (Rohwer et al., 1991; Rohwer, 1993; 
van der Werff & Richter, 1996; Li & Christophel, 
2000; Li et al, 2004). This 


therefore, be used with caution in Laureae classifica- 


character should, 
lions, although it does seem to be useful for defining 
some of the higher-level clades in the group (e.g.. the 
Litsea clade s. str.). 

The significance of other traditional characters 
often used at the generic level, such as dimerous or 
For 


example, Laurus and Neolitsea have been related 


irimerous flowers, also needs re-evaluation. 
previously based on their dimerous flowers, but not in 
our analyses. Similarly, the nature and arrangement of 
basal inflorescence involucral bracts (ie. early 
deciduous vs. tardily deciduous; imbricate vs. decus- 
sate) are traditionally important characters for delim- 
iling genera. However, whereas Actinodaphne, Para- 
sassafras, and Sinosassafras can be distinguished from 


other Laureae by the possession of early deciduous, 


imbricate bracts. these character states are conver- 
gently homoplasious in our analyses. This suggests 
that although potentially diagnostic for relating the 
latter two genera, the character is of limited 
phylogenetic value, especially given the position of 
Actinodaphne as a basal grade below Neolitsea. 
Nevertheless. the clades identified in this studv do 
provide opportunities to examine evolution of specific 
morphological characters in Laureae, a task already 


Li and Christophel (2000), and of 


particular interest will be studies of inflorescence 


initiated. by 


development. The most common inflorescence form in 
Laureae is pseudo-umbels clustered in the leaf axils, 
but our analyses show that there are apparently al 
least two different ways to produce this structure. One, 
suggested by previous studies (Li, 1985; Tsui, 1987) 
and seen in Laurus and the Litsea, Lindera, and 
Aperula | and IH clades, results from shortening of the 
internodes in the short shoot with a vegetative 
terminal bud. This vegetative terminal bud may be 
normal and scaly or naked, well developed and large, 
poorly developed, or even reduced. Sometimes the 
terminal vegetative shoot along with one or two lateral 
fertile short shoots with pseudo-umbels merge to form 
a single axillary mixed bud, and the peduncle of 
pseudo-umbels and the internodes of the short shoot 
may be developed or reduced. Given the positions of 
these taxa toward the base of the trees, this feature 
appears lo be plesiomorphic in Laureae. 

The second, derived condition, seen in the Neolit- 
sea—Actinodaphne clade, results instead from shorten- 
ing of a thyrsoid inflorescence axis and lacks a 
vegetative terminal bud. However, given the lack of 
deep-branch support for the major lineages in this 
study, definitive conclusions about the phylogeny of 
inflorescence ontogeny must await further studies. 

In conclusion. our study indicates the need for 
caution in the use of morphological similarity for 
assessing affinities between taxa in the Laureae. 
Traditional characteristics of habit, leaf venation, 
inflorescence, and floral structure appear in many 
cases to have been the result of convergent and/or 
parallel evolution and, therefore, may not be indica- 
tive of evolutionary affinity or useful for taxon 
delimitation at higher levels. Added to this is the 
possibility that some features such as pseudo-umbels 
may not be homologous. Nevertheless, the study 
identifies areas in which future research may help to 
clarify or correct problems of homology and ontoge- 
netic convergence. It provides a hypothesis for 
possible phylogenetic relationships in Laureae, albeit 
based on a single, rapidly evolving genome, and gives 
direction for future studies using multiple indepen- 
dent and possibly more conservative markers to assess 
the phylogenetic hypotheses that our results indicate. 
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This will provide the foundation for a revised 
phylogeny-based classification of the Laureae in 
which reliable synapomorphies are backed by data 
from a range of sources. 
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